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Abstract 



Recently-developed variational perturbation expansions converge exponentially fast 
for positive coupling constants. They do not, however, possess the correct left-hand 
cut in the complex coupling constant plane, implying a wrong large-order behavior of 
their Taylor expansion coefficients. We correct this deficiency and present a method 
of resumming divergent series with their proper large-order behavior. For a given set 
of expansion coefficients, knowledge of the large-order behavior considerably improves 
the quality of the approximation. 
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1) Recently, the Feynman-Kleinert variational approximation to path integrals [1] has been 
extended to a systematic variational perturbation expansion [2,3] which converges uniformly 
and fast (for the anharmonic oscillator like e - constxNl/3 m the order N of the approximation) 
[4,5]. Due to the uniformity of the convergence, this extension has given rise to an efficient 
method for extracting strong-coupling expansions from a weak-coupling expansions [6-8]. 
This has led to a new resummation procedure which makes use of a possible independent 
knowledge of expansion coefficients for weak and strong couplings [9], and should prove 
useful for extracting critical properties from lattice models of statistical mechanics, where 
both expansions are known to high order (see for example the testbooks [10]). 

In the theory of critical exponents, a different resummation procedure is required. There 
one must evaluate divergent series of which one typically knows only a few expansion coef- 
ficients plus the large-order behavior, the latter from semiclassical tunneling theory [3,11]. 
We present a method of combining the exponentially fast convergence of variational pertur- 
bation theory with the information on the large-order behavior, thus laying the grounds for 
an efficient resummation of perturbative expressions in the quantum field theory of critical 
phenomena. 

The method is completely general and holds for any physical system whose quantities, 
for instance energy eigenvalues E, possess (possibly divergent) power series expansions in 
some dimensionless coupling constant g' = g/co q {u being some energy parameter) of the 
type u' p ^2 E n (g /io q ) n , whose coefficents behave for large orders like 

E k = 1P rj+1 k rj {-a) k {pk)\ 

Such a behavior arises from a cut in the complex coupling constant plane, across which E(g) 
nas a discontinuity 

disc E( — \g\) = E( — \g\ — irj) — E( — \g\ + irj) = 2iImE( — \g\ — irj) (2) 

= (a| 5 '|/^)" (/3+1)/p e" 1/(a|fl ' l)P [l + c 1 (a\g'\) 1/p + c 2 (a\g'\) 2/p + . . .]. (3) 

A typical example is the ground state energy E(g) of the anharmonic oscillator which 
has p = 1, q = 3 and a large-order behavior 
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E k = --J-(-3/io 3 ) k T(k + 1/2)[1 - 95/72A; + ...], (4) 
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itli a discontinuity [12] 



2 ? Im£(-M " "?) = " 2 ^^y ^j e "" 3/3M[1 " ( 95 / 72 X 3 W) + ■■■]■ ( 5 ) 

In this note, we shall develop the method for this particular example and illustrate how the 
information on the large-order behavior accelerates the convergence of variational perturba- 
tion expansions. 

2) Following the procedure explained in [3], we take the weak coupling expansion of or- 
der iV, 

N 

E N = ^E n {-^)\ (6) 

n=0 

and replace u by the identical expression 

lo -> VO 2 + u 2 - fi 2 (7) 

and reexpand EJj in powers of A treating u 2 — fi 2 as a quantity of order g. The reexpanded 
series is truncating after the order n > N. 
The resulting expansion has the form 

N 

W N (g,n) = lF£a n / n (l)) (8) 

where 



71 = 1 



Forming the first and second derivatives of W^{g, fi) with respect to fi we find the positions 
of the extrema and the turning points. The smallest among these is denoted by fijv- The 
resulting Wjs;(g) = VKjv((/,0jv) constitutes the desired approximation to the energy. 
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3) The perturbation expansion of the anharmonic oscillator looks like (6) with E n = 
1/2,3/4,-21/8,333/16,-30885/128,916731/256... . The lowest-order approximation to 
the energy reads 



^.^(f + ^W^- (10) 



Extremizing this yields 



%ucosh [|acosh(flf/flfi)] g > g 

Ox = { V for (11) 
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to cos [| arccos((//(/i)] 9 < 9ii 
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ith </i = 2cu 3 £ , o/3v / 3-E'i- The result is shown in Fig. 1, where the approximation is seen to 
have a maximal error of 2% for large couplings. 

For larger negative couplings g < 0, the imaginary part of the energy is reproduced with 
the same type of error, as shown in Fig. 2. At small negative g } however, the approximation 
W\(g) has an important qualitative deficiency: It does not reproduce the imaginary part (5) 
in the interval g £ ( — </i,0) where g\ 0.2566. known from semiclassical tunneling theory. 

By going on to the approximation W^{g), the energy at positive g is found correctly 
within 0.2% (see Figs. 1 and 3). The same is true for the imaginary part at larger negative 
g (see Figs. 2, 4, and 5). Again, the imaginary part at small negative g is being missed, 
although the interval is now smaller: g £ (— #3,0) with g 0.16 

4) Let us now include the knowledge of the imaginary part at small g < in Eq. (5), 
derived via tunneling theory, into the variational perturbation expansion. First we consider 
the approximation W\(g) where the missing interval has g\ = 4:/9Sqrt3 ~ 0.2566. The 
dispersion relation for E(g) reads (with one subtraction to ensure convergence) 

where e 1 ^) starts out like 1 — (95/72)(3|g|) + . . . and looks as shown in Fig. 2. By expanding 
l/(A + g) in a power series in g } we obtain the expansion coefficients as the moment integrals 
of the imaginary part function of 1/g: 
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= J, Sa^VW5I^ /3V(9 » (13) 

We now assume only the knowledge of the leading semiclassical imaginary part (5). 
This is used to approximate the imaginary part in the entire regime where W\(g) is real for 
negative g, i.e., in the interval g £ ( — <?i, 0). There it contributes to the expansion coefficients 

These numbers are subtracted from the full expansion coefficients Ek forming E' k . For 
the new coefficients E' k , the imaginary part of W\(g) starts out at another value of 
gi = 2u 3 E' /3\/3E' 1 , for which we calculate from (13) new coefficients E' k , and so 
on. This procedure converges at g\ 0.166. For the associated coefficients E' k 
0.50117, 0.72905, -2.24059, 13.54295, -98.64571, ... we now evaluate the variational pertur- 
bation expansion W[(g). To this we add the energy associated with the dispersion integral 
(14): 

No subtraction is necessary. For positive g, the new approximation W[(g) =[ (g) + A\E(g) 
is shown in Fig. 1. It is seen to be better roughly by about 30% than the previous approxi- 
mation W\(g). 

The important qualitative advantage of the approximation Ei(g) is seen in Fig. 2. The 
imaginary part starts now at g = 0, remains constant until g = — and is approximated 
by Im W[(g) for g < —g\. 

We now go on to the approximation Wz{g) which for g > is accurate to 0.05% 
as shown in Fig. 3. The situation for the imaginary part at larger negative g is only 
slightly worse (see Figs. 4 and 5). The imaginary part is now missing in the inter- 
val g £ ( — <7i,0) with g\ —0.16. For this we evaluate the dispersion relation 
(13) and find new expansion coefficients E' k . These render a new value of etc., un- 
til the method converges at g\ 0.166. The resulting new expansion coefficients are 



E' k « 0.5000477,0.74871,-2.58993,19.84402,-214.12062,... . To the associated W^g) 
we add the contribution from the dispersion integral (14), and obtain the final result 
W^(g) =3 (g) + A 3 E(g). In Fig. 3 we see that for g > 0, the new approximation is better 
than the previous one W^{g) by roughly a factor 5. To see the convergence of the iteration 
in gi } we also plot W^(g) for the initial g\ 0.16. 

The important qualitative advantage of the new approximation is visible in Figs. 4 and 
5. There is now an imaginary part for all negative g which ensures the correct leading 
large-order behavior of the new approximation (4). The cut in the interval g £ (—0.166,0) 
is approximated by leading term in the semiclassical expression (5). 

In the special case of the anharmonic oscillator, the convergence could of course be 
improved by using our knowledge of correction factor for the imaginary part in Eq. (5). 
However, since such correction factors are generally hard to derive, we have not made use 
of them. 

5) The new resummation method incorporates into the variational perturbation expansion 
the initial, tunneling section of the left-hand cut, thus accounting for the correct large-order 
behavior of the expansion coefficients. This should prove useful in calculations of critical 
exponents of statistical systems described by quantum field theory, where the tip of the cut 
is known from tunneling theory. It will be interesting to see whether it will be possible to 
improve upon the results obtainen by Borel-Pade methods. 

Certainly, any knowledge of the strong-coupling behavior can be used to obtain better 
approximations along with the help of the procedure developed in Ref. [9]. 

Note that because of the scaling relation 



first noted by Symanzik, the resummation problem for the non-Borel-summable double-well 
potential with u 1 = — 1 can be treated in the same way [14]. 




(16) 
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1.02 




0.995 



FIG. 1. Plot of the ratio of the resummed energy E^ ppT = W[{g) (solid curve) with respect to 
the exact energy as a function of the dimensionless coupling constant g' = g/(jJ 3 for the anharmonic 
oscillator. The dashed curve shows the old approximation W\{g). The short-dashed curve indicates 
the approximation W^{g). The dotted curve indicates the approximation derived in Ref. [9] making 
use of the known strong-coupling behavior (which is not done here). 




FIG. 2. The reduced imaginary part Im W[{g) (solid curve) of the ground-state energy in 
comparison with the exact one (dotted). The curve vanishes for g £ ( — 0.26459,0), where it is 
replaced by the semi classical imaginary part. The dashed curve shows the approximation Im 
W\{g). The short-dashed curve indicates the approximation Im W^{g). 
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FIG. 3. The ratio of the approximation = W^(g) with respect to the exact ground-state 

energy E ex (g) of the anharmonic oscillator (solid curve). The short-dashed curve shows the old 
approximation W^{g), the long-dashed curve the approximation obtained for g% = 0.16 rather than 
the proper value g% = 0.166. 
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FIG. 4. The imaginary part of the approximation W^(g) (solid curve) vanishing for 
g G ( — 0.166,0) where it is replaced by the semiclassical imaginary part. The dotted curve is 
the exact imaginary part. The dashed curve is obtained for g% = 0.16. The short-dashed curve 
shows the old approximation W^{g). 
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FIG. 5. Plot of the ratio of Im W^(g) with respect to the exact Im E ex of the anharmonic 
oscillator as a function of g' = g/oj 3 (solid curve). The dashed curve shows the old approximation 
Wzig), the short-dashed curve is obtained for g% = 0.16 rather than the proper value g% = 0.166. 
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